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Abstract 

Supersonic inlet flows with mixed external-internal compressions were com- 
puted using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) me- 
thod for solving the three-dimensional unsteady, compressible Navier-Stokes equa- 
tions in conservation form. Numerical calculations were made of various flows re- 
lated to such inlet operations as the shock- wave intersections, subsonic spillage 
around the cowl lip, and inlet ‘started’ versus ‘unstarted’ conditions. Some of the 
computed results were compared with wind tunnel data. 


Introduction 

The design of the inlet system is a critical part of the design of the entire 
propulsion system. For advanced supersonic aircraft, integrated design procedures 
are required which can account for both the external and internal flow fields. Al- 
though a great deal of numerical study has been focused on purely internal flows, 
there have been few computational studies of complete inlet systems such as shown 
in Figure 1. 

The analysis of mixed compression inlet flows is complicated by the formation 
of multiple shock waves that are generated in the external flow field by the external 
compression surfaces (ramps, centerbodies, aircraft forebodies) and the cowl, and 
in the internal flow field from the cowl lip to the engine face. For subcritical flow 
(or supercritical and critical flows with flow incidence), the analysis must be able 
to treat both the external and internal flow fields to accurately predict the cowl 
spillage. Internally, the analysis must be able to treat shock wave boundary layer 
interactions and shock wave intersections. In addition, for practical inlet analyses 
the three dimensional effects due to angle of attack and nonuniform freestream must 
also be included in the analyses. 

Numerical analysis of mixed-compression inlets was led by the development of 
spatial marching algorithms based on two-dimensional method of characteristics.^ 
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Recently the characteristic method with discrete shock wave-fitting technique^ was 
combined with solutions of the boundary-layer equations to account for the vis- 
cous effects which dominate inlet flow fields.® A similar method, but combined 
with a three-dimensional swept shock interaction technique,® has been also re- 
ported. These methods demonstrated the efficient capability of computing three- 
dimensional mixed compression inlets, but were limited to the supersonic portion 
of the flow field and were incapable of treating the terminal shock and the subsonic 
flow at the engine face. 

Another approach for treating the inlet flow problem is to use the finite- 
difference solution to the full potential equation. The computational methods in- 
volve a conformal mapping technique,^ and a sequence of simpler transformations,® 
to transform the flow equation to a boundary conforming coordinate system. The 
potential equation can also be solved in Cartesian®’^® as well as in generalized^^ 
coordinates. These computational methods have been shown to be fast and stable, 
but cannot account for the important viscous effects. 

The complicated flows of mixed compression inlets are amenable to computa- 
tion by solving the Navier-Stokes equations. One of the approaches taken in this 
effort was the use of the MacCormack explicit scheme in developing a shock- 
capturing algorithm for the solution of a supersonic inviscid flow,^® and in solving 
the noniterative shock-capturing parabolized Navier-Stokes algorithm. Recently, 
the MacCormack explicit algorithm was further applied to the computations of 
two-dimensional high speed inlets*®"*^ and a three-dimensional axisymmetric inlet 
with*® and without*® angle of attack. The development of an efficient explicit code 
has been, however, hampered by the restriction on the maximum time intervals 
required for numerical stability. 

As opposed to the explicit approach, an implicit algorithm of the alternating- 
direction, implicit (ADI) finite difference scheme^®’^* permits larger time intervals 
for integration. The usage of the ADI scheme in conjunction with a diagonal implicit 
method, a finite- volume formulation^®, and implicit and hybrid methods^^ has 
resulted in the computations of a two-dimensional external compression inlet flow 
field. 

The purpose of this study is to develop an efficient numerical technique for 
accurately analyzing a three-dimensional mixed compression inlet flow field. Toward 
this end, a combined implicit-explicit approach is adopted for solving the unsteady, 
compressible Navier-Stokes equations in conservation form. The implicit part of the 
algorithm adopts the alternating-direction scheme^®’®* to integrate the continuity 
and energy equations and both the momentum equations in streamwise (^) and 
normal (t)) directions; whereas the explicit part adopts the MacCormack predictor- 
corrector scheme*^ in the azimuthal (0) momentum integration. 

The selection of the combined implicit-explicit method is based on concerns 
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over accuracy and computational efficiency of the resulting code. The inlet analysis 
requires a finely spaced grid in the boundary layer (or in the inlet wall region) to 
resolve the strong shock-wave/boundary-layer interactions and potential flow sep- 
arations. The inherently stable implicit scheme capable of accommodating large 
time steps is suitable for carrying out the ADI integration across the densely dis- 
tributed grid points in the boundary layer in the rf direction. In contrast to the 
implicit method, the explicit method (in which the numerical stability mandates the 
maximum allowable time step to be dependent on the smallest spatial scale) can- 
not efficiently handle the integration across the boundary layer. Explicit methods, 
however, can be effectively used in the azimuthal integration wherein flow variation 
is relatively minor and the fine grid clustering, as in the boundary layer, is not 
necessary. 

In the following sections, the numerical technique is described for applying this 
approach in solving the Navier-Stokes equations, the computed results are presented 
for various flow conditions, and some comparisons are made with available wind- 
tunnel data. 


Numerical Technique 


Governing Equations 

The flow is assumed to follow the full, unsteady, compressible, Navier-Stokes 
equations in conservation form: 


P,t + {p^t),t — 0 

(pK),t + ^ij,j = 0 

= 0 

where 

fl,y = P^i^j “I" P^ij 
A, = {e + p)Vi + qi- UjV) 

Here, p, V , e, and p denote density, velocity, energy, and pressure, respectively; g 
and T denote heat transfer and skin friction, respectively. 

The governing equations in cylindrical coordinates can be expressed cis 

dtV + drE + dffF + d^G = -H 
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The skin friction in terms of the cylindrical coordinates gives a form 


Trr = 2ndrUr + X(drUr + doUfflv Ur / T + dgU^) 

= 2fidzUz + \{drUr + deueir + Ur/r + dgU^) 

Tes = 2ti{dexLelr + Ur/r) -t- X{drUr + deuejr + Ur/r + dgUz) 

Trz — Tzr ~ P’i^z'^r “k ^r^z) 

Trd = Tdr ~ P-i^rUe - Ue/r + deUrjr) 

Tffz ~ z6 ~ p{d$Uz "h dgUff^ 

wherein the Stokes hypothesis A + [2/3)p = 0 is assumed. 

To obtain flow solutions on a body-conforming grid that is suitable for treating 
an arbitrary inlet geometry and implementing the necessary flow boundary condi- 
tions, a generalized transformation is applied to the governing equations using a 
form ^ = ^(z, r, d,t) , r) — r}{z, r, 6,t) , Q = &{z, r, 6,t) , and r — t . The resulting 
flow equations give the form 

+ a^Q + a©R = -H (1) 

where 

U = U/J 

P = &u + 6E + 6F + 6G 

Q = TjtU -I- TjrE -h Ttef -f rjzG 

R = 0fU + 07-C -{■ 0^F "t" 0jjG 

H = H/J 

and J denotes the Jacobian, 

+ irne^z + ienz^r - Ur)r&z - irVz^e - izVe^T 
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Here, it is denoted that E = E/J , F = F/J , and G = G/J . The first row of 
vector Equation (1) corresponds to the equation of continuity. Likewise the second 
and third rows are the momentum equations in the ^ and rj directions, respectively, 
while the fourth row is the energy equation and the fifth row is the momentum 
equation in the © direction. 

Finite Difference Procedure 

The first four rows of Equation (1) are integrated in time, subject to appropri- 
ate initial and boundary conditions, using the Beam-Warming-Steger alternating- 
direction, implicit (ADI) scheme. The last row of Equation (1), the momentum 
equation in the 0 direction, is integrated in time with the MacCormack explicit 
scheme*^ using the newly computed values (from the ADI scheme) of density, en- 
ergy, and velocities in the ^ and r/ directions. 

To apply the ADI scheme in the ^ and t) directions, the first four rows of 
Equation (1) are arranged in the following form. 


drU -|- -h d,jQt — S 


( 2 ) 


where 

S = -H, - deR, - Re~^{d(Pv + -f + H„) 

^ A. A. 

As opposed to the terms P, , Q, ,R, , and H, , which contain only the inviscid flux 
terms, the terms P„ ,Q„ ,Rv , and Ht, , contain the remaining viscous derivatives 
and heat transfer terms. 

In the Beam-Warming-Steger delta-form ADI algorithm. Equation (2) can be 
written as 

\I + hS^{A - ^^K)” - J-‘ahV|A^J) 

x|I + h6„(B - (5„L)” - A„J ](U”+^ - U'^) 

= Ar(-d^P, - + S)” - oU-i| (V^ A^)2 + (V„ A„)2 ]JU" (3) 

where second-order time differencing, h = Ar/2, and second-order central difference 
for the convection terms, 6^ and Sr,, are used in the computation. It is noted here 
that Pt, and Q„ are divided into two parts, one containing ^ and the other rj 
derivatives; that is 

P„(U,U^ ,U„) = P„i(U ,U^) + P„2(U ,U„) 

Q„(tJ,U^ ,U,) = Q„i(U,U^) + Q.2(U ,U,) 
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K and L are defined as K = , and L = Q„2 • A and result 

from the local linearization of the fluxes about the previous time level as given by 

A A A A _ _ 

A = djjPi and B = • Fourth order dissipation terms (V^A^)^ and (V, A,,)^ 

are added to the right-hand side of Equation (3) to control the numerical stability.^^ 

A A A 

The spatial derivatives and involve cross derivatives d^dr,\J^ and 

d^drfVn ; in order to use a spatially factored scheme, these terms can be time lagged 
to an n — 1 time level. These terms were neglected in the present computation, 
which reduces the order of the time-accuracy in the computation of dissipation 
terms from second to first order. For an axisymmetric body (the flow need not 
necessarily be axisymmetric), the transformations ^ = ^(2, r),»/ = 17(2, r) and 
Q = 6 give = 1 and r)e = 0 . With this transformation, the Jacobian 

is reduced to J = 0^Jz = where denotes the two-dimensional Jacobian, 
Jz = — irVz • The expressions for terms P , Q ,R , and H can also be greatly 

simplified. The left-hand-side Equation (3) can be factored in i and t) components 
and solved by sweeping in each direction by inverting tridiagonal matrices of 4 x 4 
blocks. 

The original MacCormack scheme^^ is applied to the last row of Equation (1)‘ 
or the momentum equation in © direction, using the algorithm : 

Predictor: 


n+ 1 

j,k,t 




n 



Corrector: 


=1/2(6“^., + u;X! - '•.(Py.ti - P7-U.,) - '>3(q;xi - 




where ho = At , and hi - Ar/A^, = Ar/Ar/ , and /13 = Ar/A© , and £ = jA^ , 
ri = kArf , and © = /A© . 


Grid-Generation and Boundary Condition 

Since the flow structures vary greatly in both the external and internal flow 
fields, the grid is divided into two zones; one covers mainly the external (freestream) 
region and the other the internal (inlet) region of the flow field (Figure 1). The outer 
boundary points of the external zone are specified with a uniform grid spacing. The 
inner boundary points, however, are specified with grid spacing clustered in the 
regions near the tip of the centerbody and the cowl lip, and the inlet wall in the 
throat region. 
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Once the grid is specified on its boundaries, an elliptic solver^® can be used 
to generate a smoothly spaced grid at the interior points in each zone of the flow 
field. H-grid is used in each zone to properly describe the centerbody and the sharp 
cowl lip. The interior of the grid is densely clustered in the boundary layer to 
account for the rapid flow variation and the shock-wave/boundary-layer (or shock 
wave/inlet wall) interactions. H-grids generated in the internal and external zones 
share common boundary points along a 7? = constant line extending from the cowl 
lip to the upstream boundary. 

The computations were performed on a grid consisting of two zones which 
share a strip of grid points on the common boundary extending from the cowl 
lip to the upstream boundary. The computations were carried out in each zone 
interchangeably, and the computed data in each zone were stored in a solid state 
disk for the next step iteration. Two zones of the flow field were coupled along 
the common boundary where, prior to initiating the computation on each zone, the 
boundary conditions were updated. This was done by taking the average of the flow 
quantities that had been computed on the grid points on the r} = constant lines 
directly above and below the common boundary. 

On the outer boundary of the grids, free-stream conditions are applied. Along 
the body surface of the centerbody and cowl, the tangency condition is imposed for 
the inviscid flow computations; for the viscous flow, the no-slip condition is imposed. 
An engine face static pressure at the inlet exit was imposed on the boundary using 
an explicit boundary treatment. 


Results 

An existing experimental mixed compression axisymmetric inlet model^® with 
a translating cowl was selected for analysis because of the extensive amount of data 
available for comparison. The inlet was tested in the 8- by 7-, 9- by 7-, and 11- by 
11-Foot Wind Tunnels at Ames Research Center over the Mach number range 0.6 to 
2.65. The internal contours, designed for Mach number 2.65, provided a throat area 
equal to 59 % of the capture area when the centerbody was retracted for supersonic 
operation. The experimental data were taken at a tunnel total pressure of about 1 
atm (a unit Reynolds number of about 8.53 x 10®/m at Mach number 2.65). The 
experimental data included total pressure recovery and static-pressure distributions 
along both the centerbody and cowl for various bleed configurations and mass flow 
ratios. At the design Mach number 2.65, the conical shock from the centerbody sits 
slightly ahead of the cowl lip. This Mach number was chosen for computing the inlet 
flow field at critical condition (an inlet at ‘started’ condition with terminal shock). 
In addition, Mach number 3.0 (for which no experimental data were available) was 
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chosen to study the structures of supercritical flows where the centerbody shock 
falls inside the cowl lip and the flow is supersonic throughout the inlet flow field. 
For the computation of an inlet ‘unstart’ condition, Mach number 2.5 was chosen. 

Computations and Code Performance 

The flow fields considered for computations include the flow at supercritical 
without terminal shock, critical with terminal shock, and subcritical with a de- 
tached shock. Supercritical flow without a terminal shock is supersonic throughout 
the flow field, and the computations were impulsively started with a freestream 
initial condition imposed on all grid points. Throughout the computation, zero 
pressure gradient is maintained along the downstream boundary of the external 
zone and on the inlet exit surface (engine face). Critical flow with a terminal shock 
is supersonic upstream and subsonic downstream of the terminal shock wave which 
sits just downstream of the throat. These computations were impulsively started 
with a freestream condition imposed on all grid points located outside of the inlet. 
The starting conditions for the grid points in the inlet interior were linearly varied 
from the cowl lip at the freestream supersonic condition to a subsonic condition 
of approximately M = 0.6 at the diffuser throat. From the diffuser throat to the 
engine face, the Mach number was set constant at M = 0.6 . As the computation 
was started, the static pressure at the inlet exit was raised gradually in the process 
of time iteration to a specified (experimentally measured) engine face pressure. At 
each time step, the total energy at the inlet exit was updated using the newly spec- 
ified pressure and the velocity at the inlet exit that had been extrapolated from the 
interior grid points adjacent to the inlet exit. Subcritical flow is subsonic from the 
detached shock to the diffuser throat. The computations were impulsively started 
with a subsonic {Moo = 0.8) initial condition imposed on all grid points. In the 
process of time iteration, the freestream Mach number imposed on the upstream 
outer boundaries was raised gradually to a supersonic Mach number {Moo = 2.5). 
Throughout the computations, zero pressure gradient was maintained on the down- 
stream boundary as in the case of supercritical and subcritical flow computations. 

At each iteration in time, the surface pressure was computed from the nor- 
mal momentum equation. A free stream condition was maintained on the upstream 
outer boundaries of the flow field throughout the process of integration in time. The 
computational results presented in this paper cover only an axisymmetric inlet, al- 
though the solver is capable of treating an inlet with an asymmetric configuration. 
The computations utilize a grid consisting of two zones with 138 x 16 x 49 points 
per zone in the f, 0, and rj directions (Figure 1), respectively. Inviscid flow com- 
putations on the Ames Cray XMP/48 Computer required approximately 0.24 sec 
per iteration per zone for the axisymmetric inlet at zero angle of attack, and ap- 
proximately 2.7 sec per iteration per zone for the axisymmetric inlet at an angle 
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of attack. The viscous computations took approximately 20% longer than the in- 
viscid computations. The computation takes less than 1000 iterations to attain a 
converged solution. 

Supercritical Flow Without a Terminal Shock 

Figure 2a shows the computed inviscid Mach number contours for an axisym- 
metric inlet at M^o = 3.0 and a = 0, and Figure 2b depicts the enlarged Mach 
contours of the inlet interior flow fleld. A centerbody shock wave is originated at 
the tip of the compression surface, and impinges inside the cowl lip. This shock is 
also reflected between the centerbody and cowl wall inside the inlet. From the cowl 
lip, a shock wave is originated toward the outer flow fleld. There is no shock wave 
from the cowl lip to the centerbody because the local flow angle is the same as the 
internal cowl angle. A compression wave is generated in the downstream region, 
X w 5.5, of the inlet flow field, owing to the abrupt change in the inlet geometry. 

Figure 3a shows the Mach contours of the entire inviscid flow field at Mo© = 3.0 
and a — 5°. The effects of flow incidence are shown in Figures 3b and 3c which 
depict the enlarged mach contours of the interior flow field at the windward and 
leeward planes of the inlet, respectively. The centerbody shock in the windward 
plane (Figure 3c) is closer to the centerbody than at a = 0° (Figure 2b), and 
interacts with the cowl expansion waves before impinging on and being reflected 
from the cowl surface. Two shock waves are present internally on the windward side 
of the inlet (Figure 3c); the upstream shock is the reflected shock originated from 
the centerbody, and the downstream shock is originated from the local curvature of 
the cowl inner surface. On the leeward side (Figure 3a), the centerbody shock falls 
farther upstream of the cowl lip inducing flow slippage in the lip region of the cowl. 
Internally on the leeward side of the inlet, a shock is present originating from the 
cowl lip. 

Critical Flow With a Terminal Shock 

Figure 4a shows the computed inviscid Mach contours of a critical flow (inlet 
‘started’ with a terminal shock wave) at Mo© = 2.65 and 0 = 0. Figure 4b and Fig- 
ure 4c show the computed surface pressures, along with experimentally measured 
pressure data, on the centerbody and the cowl, respectively. An experimentally 
measured^® engine face static pressure ratio, Pe/poo = 19, is imposed in the compu- 
tation as a boundary condition at the engine face. The experimental static pressure 
distributions were obtained for a bleed mass flow equal to 5.5% of the capture mass 
flow. The bleed condition was not accounted for in the computation. An abrupt rise 
in pressure, shown in both the computed and experimental results, indicates the 
presence of a shock at a location just downstream of the throat. The Mach contours 
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presented in Figure 4a show a flow which is supersonic upstream and subsonic down- 
stream of the shock. Upstream of the shock, the flow compresses supersonically up 
to the throat, and then expands from the throat to the shock. The computed throat 
shock is positioned downstream of the experimental throat-shock position owing to 
the viscous effects not being assumed in the computation. 

Subcritical Flow 

Figure 5a shows the computed Mach contours of an inviscid subcritical flow 
(inlet ‘unstarted’) at Moo = 2.5 and a = C, and Figure 5b depicts the enlarged 
Mach contours of the inlet flow field. In this case the shock is positioned upstream 
of the cowl lip, and subsonic flow spillage takes place in the region between the 
detached shock and the cowl lip. Thus, the flow condition is subsonic from the 
detached shock to the throat, an ‘unstarted’ inlet flow condition. Either an increcise 
in freestream Mach number or an increase in throat area is required to restart the 
inlet. 

Viscous Flow 

Viscous effects on inlet flow were studied by computing the supercritical flow 
case (no terminal shock) using the Baldwin-Lomax algebraic turbulence model. 
Figure 6a and Figure 6b show the results of viscous computations where the flow 
conditions. Moo = 3.0 and a = 0, are the same as for the inviscid results in Figure 2. 
Unlike the external flow, the inlet flow is greatly influenced by the boundary layer, 
whose thickness is significant relative to the inlet geometry. This is evidenced in 
Figure 6 by a large Mach number gradient across the boundary layer near the inlet 
wall. The Mach contour distribution in viscous flow (Figure 6b) is less clustered than 
that in inviscid flow (Figure 2b), indicating that the shocks are smeared because of 
the presence of the boundary layer. Also the effective inlet throat area is smaller 
and the throat Mach number is lower because of the presence of boundary layer. 


Conclusion 

A new computational method was developed that extends the scope of inlet 
aerodynamic prediction over a variety of flow conditions. 

The method was used to compute the external and internal flow field for a 
mixed compression axisymmetric inlet model. Results were obtained, at 0° angle 
of attack, for supercritical flow conditions without a terminal shock wave, critical 
condition with a terminal shock (inlet ‘started’), and subcritical condition with 
the inlet at ‘unstarted’ condition. In addition, the flow fields were computed for 
supercritical flow condition at an angle of attack of 5° . Some of these computations 
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were compared with experimental results obtained from wind tunnel tests of an 
inlet model with the same contours but those tests employed boundary layer bleed. 

The computations and experimental results compare reasonably well. The 
significant difference is in the location of the terminal shock wave; it is farther 
upstream for the experimental results because of the presence of the boundary 
layer. 

Overall, the computer code works well, but a fully viscous computational 
method is needed for the simulation of complex inlet flow fields. 
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Figure 1.- Computational grid (not all grids are shown). 
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Figure 3.— Supercritical flow without terminal shock 


effects of flow incidence, (a) Overall flow field. 
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Figure 3.- Concluded, (b) Flow field in leeward plane, (c) Flow field in windward plane. 
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Figure 4.— Critical flow with terminal shock — inlet ‘started’ condition, (a) Mach contours. 
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Figure 4.— Concluded, (b) Surface pressure on the centerbody. (c) Surface pressure on the cowl. 
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Figure 5.— Subcritical flow — inlet ‘unstarted’ condition, (a) Overall flow field, (b) Enlarged inlet flow field. 
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Figure 6.— Supercritical flow without terminal shock - viscous effects, (a) Overall flow field. 

(b) Enlarged inlet flow field. 
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